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METHOD AND APPARATUS FOR PROCESSING SIGNALS FOR RANGING 

APPLICATIONS 

10 

Field of the invention 

[0001] The present invention is related to a method 

and apparatus for receiving and processing signals, 
essentially signals of the spread spectrum type, such as 
15 are used in the Global Positioning System (GPS) . 

State of the art 

[0002] Today, satellite based positioning systems 

are widely in use. Most known are the Global Positioning 

20 System (GPS) developed by the United States Government or 
the Global Navigation Satellite System (GLONASS) developed 
by the Russian Government. The primary purpose of these 
systems is to provide position, velocity and time to any 
user on or near the Earth's surface. The user determines 

25 his position in space and time by measuring his range to at 
least 4 satellites whose position and time are accurately 
determined, the range being defined as the distance between 
the user and a satellite. 

[0003] The GPS satellites transmit signals on two 

30 carrier frequencies called LI at 1575.42 MHz and L2 at 
1227.6 MHz. The carriers are modulated by pseudo-random 

(PRN) spreading codes that are unique to each satellite, 
and by a navigation message. All satellites transmit at 
the same frequencies. The separation between the 
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satellites is possible because all the PRN codes are chosen 
to be orthogonal . 

[0004] The LI carrier is modulated by the so-called 

coarse/acquisition code ( "C/A-code" ) and the precision code 
5 ("P-code") . The C/A-code has a chipping rate of 1.023 MHz 
and a length of 1 ms, thus it contains 1023 chips (the term 
"chip" is used for the code symbols) . The same 1023 chips 
are repeated every millisecond. The P-code has a chipping 
rate of 10.23 MHz and a length of 1 week. The same pattern 
10 is thus repeated every week. The L2 carrier is modulated 
only by the P-code . The signal bandwidth at LI and L2 is 
limited in the satellite to the main lobe of the P-code 
spectrum, i.e. 20.46 MHz. 

[0005] In addition to the PRN code modulation, each 

15 of the carriers is also modulated by a 50-Hz navigation 
message, conveying all the necessary parameters to compute 
the satellite position and time. 

[0006] Signal tracking involves synchronising local 

replicas of the carrier and PRN codes of the satellite to 

20 be tracked with the incoming carrier and PRN codes from 
that satellite. This means that the receiver must have the 
ability to generate the PRN codes of all the satellites. 
Signal tracking is achieved by continuously modifying the 
phase and frequency of the carrier replica and delaying or 

25 advancing the code replica in order to maintain them locked 
on the incoming signal. Carrier tracking is accomplished 
in a phase lock loop (PLL) , while code tracking is 
performed by a delay lock loop (DLL) . By definition, when 
synchronisation is achieved, the phase of the carrier 

30 replica is the same as the incoming carrier phase, and the 
delay that had to be applied to the code replica is the 
same as that of the incoming code. 

[0007] The range measurement by a GPS receiver is 

based on the measurement of that delay, which is a direct 
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measure of the signal travel time from the satellite to the 
receiver. By multiplying the delay by the speed of light, 
the receiver computes its range to the satellite. 
[0008] Several error sources influence the accuracy 

5 of the satellite range measurement. One of the most 
cumbersome of them is the multipath error. Multipath is a 
special type of interference where the received signal is 
composed of the desired line-of -sight signal, and one or 
more constituents which have traversed slightly different 
10 paths due to reflections on surfaces or objects in the 
antenna surroundings. Multipath signals arrive at the 
receiver with a different delay, phase and power than the 
line-of- sight signal . 

[0009] The ranging error due to multipath depends on 

15 the delay, phase and power of the multipath signal with 
respect to the line-of -sight signal, and on the type of 
signal processing the receiver uses. Several digital 
processing techniques have been devised to reduce the 
effect of multipath on the ranging measurements. 

20 [0010] However, a common drawback of the available 

multipath mitigation methods is that they are not able to 
reduce multipath having a short delay. Even for the most 
recent multipath error mitigation techniques, multipath 
arriving at the receiver with a delay of less than around 

25 20 m affect the range measurement as if no mitigation 
technique were used at all. This is a serious limitation 
of the current techniques because in real life situations, 
most of the multipath signals are of short delay type. 
[0011] Document US-A-5390207 is related to a 

30 receiver for pseudorandom noise encoded signals. The 
operation of the device is based on narrowing the early- 
late spacing, to a fraction of a PRN code chip time, in 
order to mitigate multipath errors. The use of early-late 
correlators is a known technique for locking onto an 
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incoming PRN code, and the early-late spacing, i.e. the 
time spacing between early and late versions of the local 
PRN code is an important parameter when it comes to 
multipath error mitigation. The exact definition of this 
5 term is given in the next paragraph. 

[0012] Document US-A-5734674 is equally related to a 

receiver for PRN encoded signals; the receiver described in 
this document has the capability of dynamically adjusting 
the early-late spacing between correlation signals. 
10 [0013] Document US-A-5809064 is related to 

improvements to the same type of receiver as described in 
US-A-5734674, and equally based on dynamically adjusting 
early- late spacing. 

[0014] All techniques related to adjusting and 

15 essentially reducing the early-late spacing are linked to a 
modification of the DLL discriminator. A general name for 
these techniques is 'narrow spacing' , as opposed to 
techniques wherein no multipath error mitigation takes 
place ('wide spacing' techniques). It is found that prior 
20 art 'narrow spacing' techniques reduce the maximum range 
error with a factor 10 compared to 'wide spacing' 
techniques. However, an important multipath error remains 
for a wide range of multipath delays. 

[0015] Document US-A-5414729 is related to a 

25 receiver for PRN encoded signals, wherein a plurality of 
correlation signals are fed to a parameter estimator, from 
which the delay, amplitude and phase parameters of the 
direct path signal, as well as any multipath signals, may 
be estimated. This estimation however takes place by 
30 solving a system of equations, for example through a least- 
squares calculation, which requires complex and expensive 
hardware and software capabilities. 

[0016] Document US-A-5953367 is related to a 

receiver for PRN encoded signals, comprising a plurality of 



DLL correlators provided in each of the receiver's multiple 
processing channels. These correlators are combined in 
such a way as to build a DLL discriminator which is not 
affected by most of the multipath errors. This improves 
the tracking in the presence of multipath, but requires a 
significant modification of the DLL architecture, in 
comparison to wide spacing and narrow spacing DLL's. 
[0017] Document US-A-5901183 is related to a 

receiver for PRN encoded signals, wherein both the receiver 
DLL code and PLL carrier loops include a loop component 
that senses an error in its main loop caused by the 
presence of a multipath signal. A compound correlator loop 
is formed of a primary correlator and a secondary 
correlator that senses and corrects for multipath induced 
tracking error of the primary correlator. 

[0018] Document US-A-5781152 is related to a method 

wherein four pseudo random code replicas are used, namely a 
first set of early and late codes, with a first early-late 
spacing 2d, and a second set with a second early-late 
spacing 2kd. The multipath error is estimated from an 
extrapolation on the basis of correlation values measured 
with these two sets of codes. 

[0019] Document US-A-5966403 is related to a method 

and apparatus for minimizing the residual multipath 
distortion signal present in a received signal, by using a 
non-uniform weighting function for said 

estimation /minimization . 

[0020] An overview of the capabilities and 

limitations of prior art methods is given in figure 1. The 
graph shows the range error envelope for three different 
signal processing techniques. Each envelope represents the 
maximum and minimum range error as a function of the 
multipath signal delay, for a given signal-to-multipath 
power ratio; in this case this ratio is equal to 4. The 



6 

curve 100 refers to 'wide spacing' techniques, i.e. 
techniques wherein no multipath mitigation is employed. 
The curve 101 is valid for the 'narrow spacing' techniques 
referred to above, while the curve 102 corresponds to the 
5 techniques described earlier as representing the closest 
prior art. It is clear that no prior art technique is 
capable of reducing multipath errors at the lowest 
multipath delays (0-20m) . 

10 Aims of the invention 

[0021] The present invention aims to provide a 

method and a device for receiving and processing signals of 
the spread spectrum type, allowing a reduction of multipath 
errors in the lowest range of multipath delays, said 

15 reduction being in comparison to the closest prior art. 
[0022] The present invention further aims to propose 

a method and device which are technically straightforward 
and inexpensive compared to existing methods and devices. 

20 Summary of the invention 

[0023] The present invention is related to a method 

for processing at least one signal sent by a transmitter, 
said signal preferably being used for measuring the range, 
i.e. the distance between said transmitter and a receiver, 

25 said signal comprising a carrier signal modulated by a 
pseudo random noise (PRN) code, said method comprising the 
steps of : 

mixing said signal with a replica of the carrier signal, 
to acquire a baseband signal, representing said PRN 
30 code, 

multiplying said baseband signal respectively with 
M+N+l PRN code replica's (P. Mr .. f P +N ), said replica's 
being shifted in time with respect to each other, the 



value M+N+l being at least equal to four, one of said 
replica's being the punctual replica P 0/ 

calculating the M+N+l correlation values ( I_ M , I +N ) of 
said baseband signal with respect to each of said M+N+l 
5 PRN code replica's, 

calculating from said M+N+l correlation values (I_ 
m,...,I +n ), an estimate of the multipath error, said 
calculation being based on a predefined formula, 
equating said multipath error to a predefined function 

10 of said M+N+l correlation values ( I_ M , I +N ) , wherein 

said predefined formula is a linear combination of said 
M+N+l correlation values ( I_ M , I+n) * each of said values 
being normalized by the correlation value I 0 of said 
punctual replica P 0 - 

15 [0024] According to the preferred embodiment, the 

method of the invention comprises the steps of : 

mixing said signal with a replica of the carrier signal, 
to acquire a baseband signal, representing said PRN 
code, 

20 - multiplying said baseband signal respectively with three 
equally spaced replica's (P 0 , P-i, P+i) of said PRN code, 
namely an early (P-i), punctual (P 0 ) and late (P+i) 
replica, with a given early-late spacing (d) , 
multiplying said baseband signal with at least one 

25 additional replica of said PRN code, said additional 

replica being shifted in time relative to said early, 
late and punctual replica's, so that in total M+N+l code 
replica' s are used, M+N+l being at least equal to four, 
calculating the M+N+l correlation values ( I_ M , I+ N ) of 

30 said baseband signal with respect to each of the M+N+l 

PRN code replica's, 

locking the punctual code (P 0 ) to the baseband signal by 
keeping the two correlation values (I-i, I+i) between 
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said baseband signal and said early and late replica's 
(P-i, P+i) equal to each other, 

calculating the range by multiplying the delay of the 
punctual code (P 0 ) by the speed of light, 
5 - calculating from said M+N + l correlation values (I- 
M ,...,I +N ), an estimate of the multipath error, said 
calculation being based on said predefined formula, 
filtering said estimate of the multipath error and 
subtracting said estimate of the multipath error from 
10 said calculated range, yielding a corrected range value. 

[0025] In the method of the invention, said linear 

combination preferably has of the following form : 
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wherein MP represents the multipath error, d represents the 
early-late spacing, I 0 represents the correlation value of 
said punctual repli ca , I-Mf *■■/ 1+n represent the correlation 
values, oti represent M+N+l fixed values. 
20 [0026] Preferably, every one of said M+N+l replica's 

is shifted over the same time delay with respect to the 
next and/or previous replica. 

[0027] Further according to a preferred embodiment, 

said cti values are calculated according to the method 
25 comprising the steps of : 

simulating, for a fixed signal-to-multipath amplitude 
ratio and for different multipath delays : the multipath 
range error and M+N+l correlation values, 

using said simulated range errors and correlation values 
30 to obtain a system of equations, each equation equating 
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the simulated multipath range error to a linear 

combination of the M+N+l correlation values, 

obtaining said a± values by solving said system of 

equations. 

5 [0028] In this method, all except two of said a ± 

values may be set to zero. 

[0029] According to a further preferred embodiment, 

said GCi values are calculated according to the method 
comprising the steps of : 
10 - for a given signal-to-multipath SMR and for a given 
number D of multipath delays, calculating the range 
error, for multipath in phase with a simulated line-of- 
sight signal and for multipath 180° out of phase with 
said simulated line-of -sight signal, thereby obtaining a 
15 vector y comprising 2D range error values, 

calculating, for each of the 2D range errors, the M+N+l 
correlation values, and normalizing said correlation 
values by the correlation value I 0 of said punctual 
replica P 0 , to obtain a (2D X (M+N+l)) matrix C, 
20 - calculating the cti values by solving the system of 
equations : y=C.a., wherein a is a vector comprising the 
M+N+l cti values. 
[0030] The vector a may be overdetermined by said 

system of equations in which case said vector a is obtained 
25 by an optimization technique. 

[0031] According to a particular embodiment, two 

replica's are used in the estimation of the multipath 
error, and the early-late spacing d is 1/15 of a chip 
length, and the second replica P +2 is 1/15 of a chip length 
30 later than said punctual replica P 0 , and said multipath 
error estimation MP is calculated as : 



MP = -0.42 
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[0032] The invention is equally related to a 

receiver for ranging applications, said receiver comprising 
a plurality of channels for detecting and locking onto a 
5 plurality of PRN encoded signals, each channel comprising : 
a delay line, comprising M+N+l taps, M+N+l being at 
least four, for obtaining M+N+l PRN codes, one of which 
is a punctual code P 0 , one a first early code P-l, and 
one a late code P+l, with an early-late spacing d 
10 between the early and late code, 

M+N+l mixers and M+N+l accumulators to calculate M+N+l 
correlation values ( I_ M , I +N ) , 

a multipath estimator module to calculate a multipath 
error estimate MP, according to a predefined formula, 

15 said formula being a linear combination of said M+N+l 

correlation values ( I. M , I +N ) , each of said values being 
normalized by the correlation value I 0 of said punctual 
replica Po, 
a low pass filter. 

20 [0033] In a receiver according to the invention, 

said formula preferably has the form : 



M p= y a,— — 

i=-M..N ^0 — 

M 2 

wherein MP represents the multipath error, d represents the 
25 early-late spacing, I 0 represents the correlation value of 
said punctual replica , I_ M ,...,I +N represent the correlation 
values, oti represent M+N+l fixed values. 

[0034] In a receiver according to the invention, 

said oti values are preferably calculated according to the 
30 method comprising the steps of : 
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simulating, for a fixed signal-to-multipath amplitude 
ratio and for different multipath delays : the multipath 
range error and M+N+l correlation values, 

using said simulated range errors and correlation values 
5 to obtain a system of equations, each equation equating 

the simulated multipath range error to a linear 
combination of the M+N+l correlation values, 
obtaining said a± values by solving said system of 
equations . 

10 [0035] In this method, all except two of said a ± 

values may be set to zero.. 

[0036] According to a further preferred embodiment 

of the receiver of the invention, said a± values are 
calculated according to the method comprising the steps of 

15 

for a given signal-to-multipath SMR and for a given 
number D of multipath delays, calculating the range 
error, for multipath in phase with a simulated line-of- 
sight signal and for multipath 180° out of phase with 
20 said simulated line-of -sight signal, thereby obtaining a 

vector y comprising 2D range error values, 

calculating, for each of the 2D range errors, the M+N+l 
correlation values, and normalizing said correlation 
values by the correlation value I 0 of said punctual 

25 replica P 0 , to obtain a (2D X (M+N+l)) matrix C, 

calculating the a± values by solving the system of 
equations : y=C.ot., wherein a is a vector comprising the 
M+N + l cti values. 
[0037] The vector a may be overdetermined by said 

30 system of equations, in which case said vector a is 
obtained by an optimization technique. 
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[0038] According to a preferred embodiment of the 

receiver of the invention, two replica's are used in the 
estimation of the multipath error, and the early-late 
spacing d is 1/15 of a chip length, and the second replica 
P+2 is 1/15 of a chip length later than said punctual 
replica P 0 , and said multipath error estimation MP is 
calculated as : 



MP = -0.42 • 



[0039] In a receiver according to the invention, 

10 said multipath estimator module may comprise software means 
for calculating the multipath error estimate on the basis 
of a predefined formula. 

[0040] Alternatively, in a receiver according to the 

invention, said multipath estimator module may comprise 
15 hardware means for calculating the multipath error estimate 
on the basis of a predefined formula. 

[0041] The invention is equally related to a method 

for estimating a ranging error due to multipath in a 
receiver, said receiver comprising : 
20 - a delay line, comprising M+N+l taps, M+N+l being at 
least four, for obtaining M+N+l PRN codes, one of which 
is a punctual code P 0 , one an early code P-i, and one a 
late code P+i, with an early-late spacing d between the 
early and late code, 
25 - M+N+l mixers and accumulators to calculate M+N+l 
correlation values ( I-m,. • / I+n) / 

a multipath estimator module to calculate a multipath 
error estimate MP, according to the formula : 

MP= ^ a t — ^— j wherein MP represents the multipath 

i=-M..N A) ].— I/I — 



2 

30 error, d represents the early-late spacing, Io 

represents the correlation value of said punctual 
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replica , I- m ,...,I+n represent the correlation values, cti 
represent M+N+l fixed values, 
said method comprising the steps of : 

simulating, for a fixed signal-to-mult ipath amplitude 
ratio and for different multipath delays : the multipath 
range error and M+N+l correlation values, 

using said simulated range errors and correlation values 
to obtain a system of equations, each equation equating 
the simulated multipath range error to a linear 
combination of the M+N+l correlation values, 
obtaining said cti values by solving said system of 
equations . 

[0042] In this method, all except two of said cti 

values may be set to zero.. 

[0043] According to a further embodiment, said 

method comprises the steps of : 

for a given signal-to-multipath SMR and for a given 
number D of multipath delays, calculating the range 
error, for multipath in phase with a simulated line-of- 
sight signal and for multipath 180° out of phase with 
said simulated line-of -sight signal, thereby obtaining a 
vector y comprising 2D range error values, 

calculating, for each of the 2D range errors, the M+N+l 
correlation values, and normalizing said correlation 
values by the correlation value I 0 of said punctual 
replica P 0 , to obtain a (2D X (M+N+l)) matrix C, 
calculating the a± values by solving the system of 
equations : y=C.ct., wherein a is a vector comprising the 
M+N+l cti values. 

Short: description of the drawings 

[0044] Fig. 1 illustrates the ranging error 

envelopes for prior signal processing techniques. 
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[0045] Fig. 2 illustrates the signal processing in a 

receiver for a ranging application according to the prior 
art . 

[0046] Fig. 3a and 3b illustrate the effect of 

5 multipath errors on the tracking point in devices of the 
prior art. 

[0047] Fig. 4 illustrates the general concept of the 

method and device of the invention. 

[0048] Fig. 5 illustrates the optimization techique 

10 used to obtain the coefficients a±. 

[0049] Fig. 6a and 6b illustrate the principle of 

the multipath estimator module. 

[0050] Fig. 7 represents a preferred embodiment of 

the invention. 

15 [0051] Fig. 8 compares the invention to the prior 

art in terms of ranging error envelope. 

[0052] Fig. 9 represents real test data proving the 

better performance of receiver of the invention. 

20 Detailed description of the invention 

[0053] A typical implementation of the signal 

processing technique according to the prior art for ranging 
applications is represented in Fig. 2. The incoming signal 
1 is first mixed with a local replica of the carrier, 

25 produced by the carrier generator 2, in order to bring the 
signal to baseband. The frequency and phase of the local 
carrier are controlled by a PLL (not shown) . The baseband 
signal essentially represents the PRN code, modulated by 
the 50Hz navigation signal. 

30 [0054] The delay lock loop is represented by the 

group of elements 4. As known in the art, an acquisition 
step is first performed in order to acquire a 'coarse' 
match between the incoming signal and a local version of 
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the PRN code. After that, the local code needs to be 
locked onto the incoming signal. The present invention is 
related to this locking step. 

[0055] Referring again to figure 2, the local PRN 

5 code is generated in the code generator 3 at a rate which 
is continuously controlled by the DLL discriminator and 
filter 10. The code enters a delay line 5 where three 
different code replicas are generated : P 0 is the punctual 
code replica, which has to be kept aligned with the 

10 incoming code. P_i is the early replica, which is advanced 
by a fraction of a chip with respect to P 0 , and P +1 is the 
late replica, delayed by a fraction of a chip with respect 
to P 0 - The delay between any two adjacent taps in the 
delay line is the inverse of the freguency of the delay 

15 line clock 6, and is traditionally referred to as "d/2" in 
units of code chips. 

[0056] For the GPS C/A-code, the code chip duration 

is close to 1 \is, and the chip length close to 293 m. The 
delay between the early and late taps is an important 
20 design parameter, which is referred to as the early-late 
spacing, and noted y d' . Many receivers use a so-called 
"Wide spacing" of d = 1 code chip. The spacing is said to 
be narrow if d is lower than 1 chip. 

[0057] Returning to figure 2, the multipliers 7 

25 followed by the accumulators 8 compute the correlation 
between the baseband signals at the output of the carrier 
mixer 9 and each of the local code versions. The resulting 
three correlations I 0 , I-i and I+i are represented in Fig. 
3a as a function of the delay misalignment At between the 
30 incoming code and the local punctual code (P 0 ) . Fig. 3a 
corresponds to the case where no multipath signals are 
present. When Ax is positive, the local punctual code is 
late with respect to the incoming code, and the early 
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correlation is higher than the late correlation (I+i) . 

The opposite occurs when Ax is negative. As can be 
expected, the correlation between the punctual and the 
incoming code (I 0 ) reaches its peak when they are aligned, 
5 i.e. when Ax = 0 . 

[0058] The role of the DLL 4 is to keep the punctual 

code P 0 aligned with the incoming code. It is apparent 
from the figure 3a that this is achieved if the two side 
correlations I_i and I +i are equal. Therefore, the DLL is 
10 designed such that it adapts the frequency of the local 
code in such a way that the equality I-i=I+i is verified. 
As long as the DLL succeeds in doing that, the punctual 
code P 0 is locked on the incoming signal. 

[0059] This implementation gives the correct 

15 tracking point in absence of multipath. However, when 
multipath is present, the correlations represented in Fig. 
3a do not apply any more. The I 0 , I-i and I +i correlation 
values are the sum of the correlation of the local code 
with the incoming line-of -sight signal and of the 
20 correlation of the local code with the incoming multipath 
signal. The result is that the correlation profile is 
distorted. Fig. 3b illustrates the distortion of the 
correlation functions in the case of the presence of one 
multipath signal having a delay of 60 m, and a fourth of 
25 the power of the direct signal, and being in phase with the 
incoming signal. 

[0060] It can be seen that the tracking point, 

defined by the equality I +i = I_i, is no longer positioned 
at Ax = 0 . This means that the DLL is not able to 
30 correctly align the punctual code and the incoming code. 
There is a tracking error, which in Fig. 3b is around 0.1 
chips, or 30 meters. This tracking error directly 

translates into an equivalent range error. 
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[0061] The early-late spacing has a large influence 

on the sensitivity of the DLL to multipath signals. For 
example, the envelope noted "Narrow spacing" (curve 101) in 
Fig. 1 corresponds to a spacing of d = 0.1 chips. The 
5 amplitude of the horizontal parts of the error envelope 101 
in Fig. 1 is proportional to d. 

[0062] The state-of-the-art techniques to reduce the 

range errors caused by multipath rely on a modification of 
the DLL so that the deviation of the tracking point with 

10 respect to the Ax = 0 condition is kept as small as 
possible in presence of multipath. As discussed above, many 
prior art techniques are based on a reduction and/or a 
dynamic adjustment of early-late spacing d, in order to 
mitigate multipath errors. Others employ elaborate 

15 calculation techniques. 

[0063] The present invention uses a different 

approach, based on a known fact, namely that the signal 
amplitude measurement (or equivalently the carrier-to-noise 
ratio measurement) reported by a receiver is highly 

20 correlated with the range error caused by multipath. This 
is disclosed in the document "Multipath Mitigation, 
Benefits from Using the Signal-to-Noise Ratio", J.M. 
Sleewaegen, Proceedings of the ION GPS-97 Meeting, pp. 531- 
540, 1997. One of the key properties of the signal 

25 amplitude measurement is that it is most sensitive to short 
multipath delays. The present invention makes use of this 
property to derive a multipath error estimator that is 
operating even for short multipath delays. 

[0064] In the invention, a conventional narrow 

30 spacing DLL is used for the tracking, i.e. a DLL with a low 
early-late spacing, for example d = 1/15 of a chip length. 
Tracking is thus performed in the classic way, i.e. by 
keeping the early and late correlation values equal to each 
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other (see figure 3b) . At least one additional correlation 
value is calculated, based on at least one additional 
replica of the PRN code, shifted in time with respect to 
the punctual, early and late versions. The multipath 
5 errors are estimated a posteriori by an independent 
multipath estimator module, on the basis of a predefined 
formula, comprising at least the correlation value between 
the incoming signal and said at least one additional PRN 
code replica. The method of the invention yields a minimum 
10 impact on the tracking process, and allows to easily turn 
on or off the multipath estimation process, without 
modifying the tracking process. 

[0065] Figure 4 illustrates the method of the 

invention, and shows equally the characteristic building 

15 blocks that need to be present in a receiver according to 
the invention. In the conventional implementation shown in 
figure 2, the correlation peak is measured at three 
different points, by providing an early, late and punctual 
version of the code. In the proposed invention, the 

20 correlation peak is measured with respect to the same three 
versions of the code, and in addition to at least one more 
version, shifted in time with respect to the first three. 
[0066] According to the general case presented in 

figure 4, M+N + l correlation values are computed, namely M 

25 early versions ( P_ M , P_x) , one punctual version (P 0 ) , and N 
late versions ( P +1 , P +N ) . This is done by using a delay 
line 20 having M+N+l taps, commanded by a delay line clock 
21 and generating M+N+l versions of the local code. The 
delay between the taps is still noted d/2, and the delay 

30 between tap i and tap 0 is id/2 ( i = -M, +N) . 

[0067] As in the prior art receiver, a carrier 

generator 2 and code generator 3 are used, and a plurality 
of mixers 7 and accumulators 8. All these elements in 
themselves are identical to the ones used in the prior art 
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receivers. For the tracking, a conventional DLL 

discriminator and filter 10 is used, to perform tracking in 
the classic way, namely based on one early (P-i) , one 
punctual (P 0 ) and one late PRN version (P+i) . 
5 Independently, all the M+N + l correlation values ( I. M , I +N ) 
are fed into the multipath estimator module 23 at each 
calculation of the range. This module derives M+N+l 
independent estimates of the signal amplitude by scaling 
each of the M+N+l correlation values by 1/ ( 1- | i | d/2 ) , 

10 i=(-M, ..,N). When multipath is present, each of these 
signal amplitude estimates exhibits an error which is 
highly correlated with the range error, although it is 
different for each estimate. The invention lies in taking 
advantage of these differences to build a multipath range 

15 error estimator. More specifically, it will be shown that 
an appropriate linear combination of the M+N+l signal 
amplitude estimates, normalized by the punctual correlation 
value I 0 closely matches the ranging error due to 
multipath. In other words, the range error due to multipath 

20 may be closely approximated by the following formula: 



MP 



NT 1 /, 

i=-M..N J Q — 



This estimation of the multipath error MP is characteristic 
25 to the invention. The coefficients a± are constant. They 

are computed only once, during the design of the receiver. 

After the calculation according to the predefined formula, 

the noise on the estimation (MP) is filtered out by a low 

pass filter 24. The result is then subtracted from the 
30 range derived from the early-late tracking by the DLL, 

resulting in largely removing the multipath error from the 

range measurement. 
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[0068] This predefined formula allows a very fast 

estimation of the multipath error, compared to existing 
techniques, in particular compared to the technique 
presented in US5414729 wherein the multipath error is found 
5 at each ranging step, by solving a system of equations. In 
the method of the invention, no estimation is made of 
multipath parameters such as delay, phase and amplitude. 
Only the multipath error itself is estimated. 
[0069] The coefficients oci are obtained by a 

10 calculation described in hereafter. Multipath range errors 
and M+N + l correlation values are simulated for a predefined 
signal-to-multipath amplitude ratio (SMR) , and for 
different values of multipath delays. Formulas to perform 
this computation are well documented, for instance in the 

15 document "Multipath Mitigation, Benefits from Using the 
Signal-to-Noise Ratio", J.M. Sleewaegen, Proceedings of the 
ION GPS-97 Meeting, pp. 531-540, 1997. This simulation can 
be easily performed using a known calculation software such 
as Matlab. 

20 [0070] The a ± f s are found by equating the 'known' 

simulated multipath range error (expressed in code chips) 
to a linear combination of the normalized correlation 
values. The coefficients a ± then become unknowns in a 
system of equations. It is also possible to set one or 

25 more cci' s to zero, thereby excluding the corresponding 
correlation (s) from the linear combination. At least two 
oti values however must not be put to zero. 

[0071] According to the preferred embodiment, the a± 

coefficients are computed using the following method: 
30 The multipath error profile is simulated for one particular 
value of the signal-to-multipath amplitude ratio (SMR) , for 
example 10, and for the worst case conditions where 
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multipath is in phase with the direct signal, and 180 
degrees out of phase. 

[0072] To cover the whole range of multipath delays, 

it is sufficient to sample the multipath profile at a given 
5 number of D points linearly equally spaced at multipath 
delays between for example 0 and 1.2 chips. When D=50, 
this yields 100 samples of the multipath profile: 50 for 
the in-phase case, and 50 for the out-of-phase case. This 
is shown in figure 5 : curve 25 is defined by the samples 
10 26. 

[0073] Let y be the vector of 100. multipath error 

samples. For all the corresponding 100 multipath 

conditions, the correlation values at taps i (i = -M..N) 
are computed, and normalized by the factor 1/I 0 / ( 1- I i I d/2 ) . 
15 This yields M+N+l vectors of 100 normalized correlation 
values, noted Cj ( j = l . . M+N+l ) . Let C = [c 2 c 2 . . . c M + ml ] be-, 
the matrix [100 X (M+N+l)] of correlation samples. 
[0074] The a ± coefficients are chosen to best fit 

the vector of multipath estimates MP = J]a .c y to the vector 

j=-M ..N 

20 of multipath samples y. 

In other words, the coefficients oti are the solution of the 
system Ccc=y. As a is overdetermined by this set of 
equations, given that N+M+l is normally smaller than 100, 
the solution is obtained by an optimization technique, 

25 preferably a least-squares optimization. By removing one 
or more columns from the C matrix, it is possible to 
exclude one or more correlations from the linear 
combination, or equivalently to set the corresponding oti 
coefficient to zero. 

30 [0075] The curve 27 in the figure 5 represents an 

example of the vector MP after least-squares fitting to the 
multipath samples of curve 25. The residual error is shown 
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as curve 28. This is the remaining error after application 
of the invention. In the above description, the values of 
M and N and the dimension of the vector y can be chosen. 
The larger these values, the more accurate will be the 
5 curve fit. 

[0076] The optimisation of the cci coefficients has 

been performed for a particular value of the SMR, e.g. 10. 
The invention lies in the finding that computing oci for 
other values of SMR yields very similar values. In other 
10 words, the a± coefficients differ little for a wide-range 
of SMR. This is essential to the invention as it is 
impossible to predict the SMR without concrete knowledge of 
the reflecting obstacles. 

[0077] For instance, taking the example of M=0, N=3 

15 and d=l/15, the a coefficients are shown in the following 
table for different values of the SMR ranging from 10 to 
100. It is assumed that this range of SMR is the most 
prevailing in a real-life multipath environment: 



SMR 


a 0 


a. 


0.2 


a 3 


100.0000 


0. 0207 


-0 . 6541 


0. 6618 


-0. 0283 


77 . 4264 


0. 0196 


-0. 6523 


0. 6606 


-0. 0280 


59. 9484 


0. 0182 


-0. 6498 


0. 6591 


-0. 0275 


46. 4159 


0. 0164 


-0. 6466 


0. 6571 


-0.0269 


35. 9381 


0 . 0140 


-0 . 6423 


0. 6544 


-0. 0261 


27 . 8256 


0 . 0110 


-0. 6367 


0. 6507 


-0. 0249 


21 . 5443 


0. 0069 


-0. 6292 


0. 6456 


-0. 0233 


16. 6810 


0. 0017 


-0. 6191 


0. 6386 


-0.0211 


12 . 9155 


-0. 0054 


-0. 6042 


0.6264 


-0.0168 


10.0000 


-0. 0147 


-0.5836 


0. 6089 


-0.0104 
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The table shows that the absolute variation of the cti's 
over the selected SMR range is small. 
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[0078] To be specific, the following linear 

combination is obtained for M=0, N=3, d=l/15, and based on 
a SMR of 10: 



MP = -0.015 — - 0.584— 1 + 0.609 7 " 2 1 -0.01^- 1 



/„ I 0 \-dl2 I 0 \-d I 0 \-3d/2 

The case of M=0 is equivalent with the situation wherein 
M=l and a. x is set to zero. In a receiver of the 
invention, there will always be at least one early 

10 correlator, one late correlator and one punctual correlator 
(and at least one additional correlator) , so that M and N 
are at least equal to one, when indicating the number of 
physical early and late correlators in the receiver. The 
fact that one or more of the a ± values may be set to zero 

15 however allows in the calculation of the range error, to 
disregard one or more correlation values, in this case I_i, 
so that M=0 for the purpose of calculating MP. 
[0079] Fig. 6a and 6b compare the ranging error due 

to multipath and the estimation of it from the above 

20 formula, as a function of the multipath delay, for two 
different multipath signal amplitudes: 10 and 100 
respectively. In figures 6a and 6b, the upper curve 
corresponds to a multipath component arriving in-phase with 
the line-of-sight component, the lower curve correspond to 

25 a 180° phase shift. It is apparent that the multipath 
estimate (curves 30) given by the linear combination given 
above for M=0 and N=3, closely matches the range error 
(curves 31), even for very short multipath delays. The 
curves 32 show the residual .range error. The cti 

30 coefficients in the above formula have been optimised to 
best approximate the error from multipath having an SMR of 
10 (case of Fig. 6a,). However, the same set of 
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coefficients provides a pretty good approximation of the 
error for other multipath amplitudes (see Fig. 6b, signal- 
to-multipath amplitude ratio = 100) . 

[0080] We can therefore summarize the inventive 

effect of the invention as follows : 

the calculation of the a ± coefficients is insensitive to 
the SMR value chosen for this calculation. In other 
words, whether the calculation is based on a SMR of 10 
or 100, the cti values do not change much. 
- With a given set of a a ± values, calculated on the basis 
of a particular SMR (e.g. 10), it is possible to 
estimate multipath errors of a wide range of multipath 
signals, i.e. signals having a wide range of the SMR. 
This allows to conclude that the a ± values can be 
calculated once during the design of the receiver, based on 
one SMR value, within the realistic range. An acceptable 
estimation of range errors due to multipath of different 
SMR, will be possible using this ONE fixed formula. 
[0081] The invention is equally related to a 

receiver, comprising a multipath estimator which calculates 
the range error based on the linear combination with fixed 
cti values as disclosed above. 

[0082] The design parameters that have to be 

optimised in the design of the receiver are the number of 
correlators and a ± coefficients, M+N+l. They result from a 
trade-off between receiver complexity, accuracy of the 
multipath error estimation, and noise on the estimation. 
Generally speaking, using more correlators (needing to 
increase the number of taps of the delay line) results in 
being able to better estimate the multipath error, at the 
expense of receiver complexity, and noise. 

As a result, a receiver according to the invention 
comprises a plurality of channels, each channel being able 
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to detect and lock onto a different PRN encoded signal 
(from a different satellite) . Each channel comprises the 
elements shown in figure 4 : a carrier generator 2, carrier 
signal mixer 9, delay line 20 and delay line clock 21, 
local code generator 3, the mixers 7 and accumulators 8 for 
all PRN versions P-„, P +M , as well as a DLL discriminator 
and filter 10, designed for "narrow spacing' type tracking 
(d < 1 chip) . The discriminator and filter may be software 
based applications, producing command signals to be used as 
inputs for the local PRN code generator 3. A receiver 
according to the invention is characterised by the presence 
of the multipath estimator module 23 and the low pass 
filter 24, wherein the multipath estimator is a means for 
estimating the multipath ranging error on the basis of a 
linear combination with predefined fixed cci coefficients. 
It is to be understood that these coefficients are 
preferably obtained during the design of the receiver by 
the calculation described above. However, small deviations 
from the optimized values of a± will still yield good 
results. This multipath estimator module can be a software 
application, that performs the action of calculating the 
estimated multipath error on the basis of the correlation 
values (I- m ,...,I+n) and the above formula. The estimator may 
equally be performed in hardware, such as a semiconductor 
chip. Also the filter 24 may be a software application or a 
hardware application. 

[0083] The invention is equally related to the 

method per se, already described, of calculating the a± - 
values on the basis of a simulated signal. 

Description of a preferred embodiment of the invention 
[0084] Simulations have shown that a good compromise 

between accuracy of the multipath error estimation and 



10 
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noise on the estimation can be achieved by using only one 
additional correlator, at position +2 in the delay line 
(M=l, N=2), i.e. based on a version of the PRN code that is 
delayed over a time equal to d with respect to the punctual 
version. The preferred design is represented in Fig. 6. 
The components printed in bold are the additional 
components with respect to the conventional receiver. 
The multipath estimator 23 of figure 7 computes the 
following function of its inputs: 



f i i \ 



MP = -0.42 



1- 



1+2 1 
I 0 1-d 



The coefficient -0.42 is the result of the optimisation of 
the oti parameters described above so that the resulting 

15 multipath estimation best fits the actual multipath error 
in the least-squares sense, for a SMR of 10. In the 
preferred design, the parameter d is set to 1/15. The 
values cc-i and oc+i have been set to zero in the calculation 
to obtain the above formula. 

20 [0085] The estimation is then fed to a low pass 

filter 24. In the preferred design, the low pass filter 
has a noise equivalent bandwidth B n of 0.1 to 1 Hz. 
[0086] Fig. 8 presents a comparison between the 

multipath envelopes obtained using a conventional DLL with 

25 a narrow spacing of d=l/15 chips (curve 101), the 
techniques representing the closest prior art (curve 102) 
and the proposed new technique (curve 103) . The curves are 
relative to a signal-to-multipath amplitude ratio of 10. 
It can be seen that the proposed technique yields the best 

30 results for multipath delays shorter than about 20 m. 
Although the existing techniques perform better for medium 
to large multipath delays, the new technique yields better 
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results in real life situations because most of the 
multipath signals fall in the short delay region. 
[0087] Also, contrary to prior art techniques, the 

multipath estimation is filtered independently from the 
5 range measurement in the low-pass filter 24, allowing to 
use a low noise equivalent bandwidth B n , and hence to keep 
the noise on the estimate low. This makes sense because 
multipath errors typically only contain very low frequency 
components. As an illustration, it can be demonstrated that 
10 the standard deviation of the noise on the preferred 
multipath error estimator is given by: 



N 0 d(2-d) 



" S \ (l-d) 2 

where X c = 293m for the GPS C/A-code. As an example, <j M p is 
only 0.1m for d=l/15, B n =0.1Hz and a nominal C/N 0 of 45dB- 
15 Hz. 

[0088] Finally, some real measurement results are 

shown in figure 9, which is showing pieces of real 
multipath signals observed during a one-day long 
experiment. The curve 40 shows the error affecting the 

20 narrow correlator tracking, and the estimation through the 
method of the invention, is shown by curve 41. The 
corrected measurement is shown by curve 42. Curve 43 shows 
the corrected measurement from a receiver using the strobe 
correlator technique disclosed in US-A-593367. Curve 43 is 

25 drawn underneath curve 42, but should actually be 
superimposed onto curve 42. The present way of drawing the 
curves allows to compare the shape of curves 42 and 43. 
From this comparison, it is apparent that the curve 42 
(invention method) deviates much less from the horizontal 

30 than curve 43. This means that the average range error is 
smaller using the invention method. Figure 9 illustrates a 
typical case where long-period multipath occurs, which are 
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often associated with short delays. This is exactly the 
case wherein the present invention delivers improved 
results . 

[0089] A complete analysis of the 24-hour data 

5 demonstrated that the method of the invention achieved a 
reduction of the multipath error by 20% with respect to the 
narrow-correlator technique, and 5% with respect to the 
strobe correlator . 



